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Abstract. Scattering differences induced by frozen particle microphysical properties are investigated, using the vertically (V) 
and horizontally (H) polarized radiances from the Global Precipitation Measurement (GPM) Microwave Imager (GMI) 89 and 
166 GHz channels. It is the first study on frozen particle microphysical properties on a global scale that uses the dual-frequency 


microwave polarimetric signals. 


From the ice cloud scenes identified by the 183.3+3 GHz channel brightness temperature (TB), we find that the scattering by 
frozen particles is highly polarized with V-H polarimetric differences (PD) being positive throughout the tropics and the winter 
hemisphere mid-latitude jet regions, including PDs from the GMI 89 and 166 GHz TBs, as well as the PD at 640 GHz from the 
ER-2 Compact Scanning Submillimeter-wave Imaging Radiometer (CoSSIR) during the TC4 campaign. Large polarization 
dominantly occurs mostly near convective outflow region (i.e., anvils or stratiform precipitation), while the polarization signal is 
small inside deep convective cores as well as at the remote cirrus region. Neglecting the polarimetric signal would easily result 
in as large as 30% error in ice water path retrievals. There is a universal “bell-curve” in the PD - TB relationship, where the PD 
amplitude peaks at ~ 10 K for all three channels in the tropics and increases slightly with latitude (2-4 K). Moreover, the 166 
GHz PD tends to increase in the case where a melting layer is beneath the frozen particles aloft in the atmosphere, while 89 GHz 
PD is less sensitive than 166 GHz to the melting layer. This property creates a unique PD feature for the identification of the 
melting layer and stratiform rain with passive sensors. 

Horizontally oriented non-spherical frozen particles are thought to produce the observed PD because of different ice 
scattering properties in the V and H polarizations. On the other hand, turbulent mixing within deep convective cores inevitably 
promotes the random orientation of these particles, a mechanism works effectively on reducing the PD. The current GMI 


polarimetric measurements themselves cannot fully disentangle the possible mechanisms. 


1. Introduction 


Cloud processes play an instrumental role in determining the surface precipitation characteristics (Rutledge and Hobbs, 1983; 
Gedzelman and Arnold, 1994; Igel et al., 2013; Tao et al., 2013). In particular, cloud ice processes are arguably among the most 
poorly understood, in part due to various ice crystal types and sizes, as well as multiple pathways in ice particle formation and 
evolution. While ice microphysical processes themselves cannot be readily observed from space, the integrated effects of these 
processes (e.g., cloud and precipitation structures, microphysical/macrophysical properties) can be inferred using remote sensing 
techniques (e.g., passive, active, multiple instruments, etc.). 

It is imperative to distinguish between ice cloud and frozen precipitation in the atmosphere although they are both composed 
of various non-spherical particles (Wallace and Hobbs, 2006). Weather and climate models treat these hydrometeors quite 
differently by suspending ice clouds for a long duration but removing frozen precipitation unrealistically fast (Li et al., 2013). 


The models account ice clouds in radiation calculations but often ignore the radiative fluxes and heating rates from frozen 
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precipitation in the air (Waliser et al., 2011). Hence, remote sensing of microphysical properties in frozen precipitation and their 
connection to ice clouds above and surface precipitation, like this work, will provide a valuable surrogate on climate models in 
representing cloud-precipitation processes. 

Retrieving microphysical properties of frozen precipitation has been a challenge for spaceborne remote sensing. It depends not 
only on macrophysical variables (e.g., column integrated mass amount, particle size distribution, etc.), but also on frozen particle 
shape and orientation (Comstock et al., 2007). Small pristine ice particles are observed with no preferred orientation when the 
ambient temperature is very cold, but become more horizontally oriented in warmer ice clouds (Noel and Chepfer, 2010). On the 
other hand, for large ice particles such as snow aggregates, it is often difficult to determine their orientation due to irregular 
shapes that are subject more to aerodynamic conditions than temperature (Xie, 2012). 

Passive and active microwave techniques at high-frequencies (> 89 GHz) are proven to be quite valuable for ice cloud and 
snowfall remote sensing because of dominance of ice scattering, whereas low frequency microwave channels are more suitable 
for detecting liquid precipitation and water clouds (Skofronick-Jackson and Wang, 2000; Wu and Jiang, 2004). To infer the 
particle shape and orientation properties, polarimetric measurements are needed and have been explored in a number of earlier 
studies. For example, Czekala [1998] pointed out the possibility of using off-nadir paired polarized channels at 200-340 GHz to 
measure horizontally oriented ice crystals. Recently, Miao et al. [2003], Xie and Miao [2011] and Xie et al. [2012] applied a 
radiative transfer model to ground-based measurements, and found up to -10K polarization differences between the vertically (V) 
and horizontally (H) polarimetric (V-pol and H-pol, hereafter) observations at 150 GHz during snowfall events. For spaceborne 
remote sensing, Prigent et al. [2001; 2005] analyzed the polarimetric observations at 37 and 85 GHz from SSMI/S (Special 
Sensor Microwave Imager Sounder) and 85 GHz from TRMM-TMI (Tropical Rainfall Measurement Mission Microwave 
Imager), where they attributed significant polarimetric differences to horizontally oriented non-spherical liquid or frozen 
precipitating particles (Prigent et al., 2001; 2005). Defer et al. [2014] found similar amplitudes of polarimetric signal at 89 GHz 
and 157 GHz using the MADRAS (Microwave Analysis and Detection of Rain and Atmospheric Structures) instrument onboard 
the Megha-Tropiques Mission. Studies using ground-based radar also suggested similar polarimetric difference distribution (e.g., 
Homeyer and Kumjian, 2015). Davis et al. (2005) studied the 122 GHz polarimetric radiances at a very large oblique view angle 
and found noticeable differences between V- and H-pol measurements, supporting the global presence of horizontally-aligned ice 
particles in the upper troposphere. In a further study with Monte-Carlo simulations, Davis et al. (2007) showed such a 
polarimetric difference remains detectable at 190-664 GHz from large oblique angles. In addition to the microwave techniques, 
lidar polarimetry at visible channels has also been used to infer ice particle orientation (e.g., Hu, 2007; Noel and Chepfer, 2010; 
Zhou et al., 2012, 2013), but limited to the top layer of ice clouds due to poor cloud penetration with lidar techniques. 
Microwave polarimetric observations at 90-200 GHz are capable of penetrating thick clouds, and yet very much under-utilized. 
With careful consideration of surface or liquid cloud effects, these polarimetric measurements can provide useful information on 
the shape and orientation of frozen particles in the atmosphere. 

The launch of Global Precipitation Measurement (GPM) core satellite in February 2014 enables new investigations on 
microphysical properties of frozen ice particles and their connection to surface precipitation. Conically polarimetric scanning at 
high microwave frequencies (89 and 166 GHz) from GPM’s Microwave Imager (GMI) enhances sensitivity to frozen 
precipitation particle scattering and preserves the polarimetric information from ice clouds. Combining with GMI radiances at 
low microwave frequencies and GPM radar reflectivity, simultaneous retrievals of frozen particle microphysical properties (e.g., 
bulk size and orientation parameters) can be achieved on a global basis. 

In this paper we present an analysis of the global GMI’s polarimetric measurements at 89 and 166 GHz to infer the 


microphysical properties of frozen particles above precipitating systems. In the analysis we also include the airborne 640 GHz 
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polarimetric observations from the NASA’s Compact Scanning Submillimeter-wave Imaging Radiometer (CoSSIR) instrument 
to extend the sensitivity to small ice crystals. Horizontally oriented non-spherical frozen particles, and the mixing induced by 
vertical motion are thought as the leading causes of the observed radiance polarimetric differences. Radiative transfer models 
(RTMs) are employed to further quantify these hypotheses. 

This paper is organized to describe the data analysis technique and radiative transfer models in Section 2, followed by the 
GMI and CoSSIR observation results in Section 3. The working hypotheses on the observed polarized signals are discussed in 


the subsequent section, with conclusions and future directions in the end. 


2. GMI data, analysis, and radiative transfer models 
2.1 GMI polarimetric measurements 


The GPM core satellite consists of a Dual-frequency Precipitation Radar (DPR) and a passive GMI instrument with 13 
channels between 10 to 190 GHz, among which the 10.65, 18.7, 36.5, 89 and 166 GHz channels have V- and H-polarizations. In 
this study, we use only the 89 and 166 GHz polarization data because at these high frequencies the scattering contributions from 
ice particles and frozen precipitation become significant. The GPM core satellite flies at an altitude of 407 km in a non-Sun- 
synchronized orbit, covering a latitude range of 68°S to 68°N. Its slow progressing rate over local time makes it feasible to study 
diurnal variability of cloud and precipitation. Its wider-than-TRMM latitude coverage now allows investigations of 
cloud/precipitation properties in the extra-tropics. GMI has a forward conical scan off-nadir at an angle of 48.5° (52.8° incidence 
angle at the surface). GMI contains two sets of footprint sizes with different scan swaths. For low-frequency (LF) channels 
below and including 89 GHz, the scan swath is 885 km, but it becomes 835 km for high-frequency (HF) channels. Such a 
disparity in footprint size and swath requires a post-processing and regridding the raw radiance measurements (L1B) in order to 
intercompare the 89 and 166 GHz polarimetric measurements. Thus, we choose to use L1C-R product instead of L1B in our 
study. More details on comparing the 89 and 166 GHz datasets can be found in the Appendix A. 

In this study we define the polarimetric radiance difference (PD) as ATB = TBy-TBy. Fig. | shows a squall line event 
revealed in the GMI’s dual-pol 89 and 166 GHz radiances and PDs on April 29, 2014. The squall line exhibits a clear centerline 
of deep convective cells and periphery of anvil clouds in the radiance maps (Fig. 1a and 1b), where brightness temperature (TB) 
is depressed strongly in both channels but their PDs are relatively small along the deep convective line for both channels. Away 
from this convective center line, PDs maximize in the anvil/stratiform precipitation region (on the order of 10 K) before diminish 
gradually in the remote clear-sky/cirrus regions (0-5 K). This variation from convective cores to anvil outflow and further to 
clear sky is more obviously seen in the 166 GHz than 89 GHz maps, because of the increasing contribution of ice particle 
scattering at higher frequency microwave channels. It is also evident in Fig. 1c, where the 89 GHz PDs are more sensitive to 
polarimetric signals from ocean surfaces than the 166 GHz, exhibiting a large PD contrast between land and ocean (notice that 
the dark-red colored ocean observations in Fig. lc actually way exceeded the colorbar range. In order to emphasize the common 
feature and disparities of the PD distribution, we chose to keep the same colorbar for Fig. 1c and 1d though). Such a land-ocean 


contrast is not readily seen in the TB maps. 
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Figure 1: A squall line event on April 29, 2014 passing through the Southeast U.S. as observed with GMI 89 GHz (left panels) 
and 166 GHz (right panels) polarimetric channels. The GMI observed radiance distributions and polarimetric differences (PDs) 
are shown in the top and bottom, respectively. The dark-red area over the ocean in panel (c) greatly exceeds the color scale with 
values mostly in the range of 16 — 29 K. They are presented by the black open dots in Fig. 2. 

Capable of penetrating clouds though, the PDs from 89 and 166 GHz contain surface polarization signals that need to be 
taken into account prior to cloud analyses. The polarized surface contributions are particularly strong over water where the 
degree of linear polarization (DoLP) of microwave radiances depends on surface roughness induced by wind and sea foams 
(Meissner and Wentz, 2012, Chen, 2014). As seen in Fig. lc, the 89 GHz PDs exhibit a large contrast between land and ocean, 
but such a contrast disappears at 166 GHz. This feature is again manifested in Fig. 2, where observations over ocean, marked by 
black open circles, spread to large PDs at the warmest end from 89 GHz observations, while this branch is not distinguishable at 


166 GHz. 
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To study scattering properties of cloud and frozen particles jointly with the 89 and 166 GHz, we need to isolate cloud cases 


from the clear sky where the surface contributions are not negligible. To flag heavy cloudy cases, we use the simultaneous 


183.343 GHz measurements, which show little sensitivity to surface emission due to the strong water vapor absorption but still 
with enough sensitivity to deep convective clouds. As described in Gong and Wu (2013), a “30 method” can be used to 
distinguish between clear- and cloudy-sky scenes for the Microwave Humidity Sounder (MHS) tropical observations. The “3o 
method” first identifies the peak (TBo) and standard deviation (co) of the Probability Density Function (PDF) of a month-long 
radiance observation in a given region or a latitude belt, and then the cloudy-sky scene is defined at places where TB < TBy - 30. 
In this study, we apply the same methodology to select cloudy-sky scenes. Cloud occurring frequency after applying this “30 
method” is comparable to that observed by CloudSat radar (see Appendix D for a detailed map). Because of ocean/land 
differences, separate PDFs are computed for ocean and land surfaces in each 10° latitude bin. The three GMI channels (89, 166 
and 183 GHz) have roughly the same (7.2 km x 4.4 km) footprint size, which minimizes the sampling error due to beam-filling 
effects between different frequencies. The combined sensitivities of the three channels will help to distinguish the scattering 
differences between frozen precipitation (large particles) and ice clouds (small ice crystals). 

As a caveat, GMI swaths are mismatched between the LF and HF channels, offset by as much as 4 degrees in pointing (LIC, 
ATBD). Therefore, the L1C data are not appropriate for joint analyses of 89, 166 and 183 GHz measurements. Instead, we use 
the LIC-R product that is derived from the L1B radiances using the nearest-neighborhood matching method to register 
collocated LF and HF observations. This re-registration procedure eliminates ~ 7% of the data, of which most occur at the edge 
of the LF scan swath. These missing data do not impact our climatological study. A comparison of the L1C-R and L1B datasets 
is given in the Appendix A. 

While the 89 and 166 GHz polarimetric measurements are useful to infer microphysical properties of frozen precipitation, PD 
at a higher frequency can reveal cloud properties with smaller ice crystals. To expand our study to ice clouds with small ice 
crystals (e.g., cirrus), we also include an analysis of the 640 GHz V- and H-pol data acquired by NASA airborne Compact 
Scanning Submillimeter-wave Imaging Radiometer (CoSSIR) during the Tropical Composition, Cloud and Climate Coupling 
(TC4) campaign in July-August 2007 (Evans et al., 2005; 2012). The CoSSIR instrument had a similar conical scan geometry to 
GMI in both forward and backward directions. About ~39000 CoSSIR conical samples from three TC4 flights are analyzed. 
These samples were mostly over the Pacific Ocean near Central America and large enough to commensurate with the GMI 
statistics over the tropical oceans. We found that CoSSIR measurements were unreliable during the ascent/descent flight in 


altitude, and therefore only measurements at a cruise altitude above 19 km are used. 


2:2 Data analysis 


PD-TB, relationship is a very useful diagnostic tool and will be used throughout this study. Fig. 2 shows a scatter plot of 
PD vs. TB, from the same squall line event reported in Fig. 1. As expected, for large oceanic PDs seen at 89 GHz, there is a 
distinct branch in Fig. 2a, showing a strong linear PD-TB, relationship at the warm or clear-sky TBs that corresponds to the 
surface polarization signals. This branch is separated from the frozen particle scattering cluster that exhibits a bell curve in the 
PD-TB, relationship. The bell curve is characterized by a peak of PD = ~10 K at TB, = 220 K. Note that the surface branch is 
absent at 166 GHz due to strong water vapor absorption in the subtropics (Fig. 2b), but the bell curve is similar to the one at 89 
GHz with peak of PD = ~10 K at TB, = 190 K. The common bell curve indicates that PD is small at the warmest (i.e., thin cloud 
or clear-sky) and coldest TB (i.e., deep convection) cases. As shown in the later sections, the bell curve appears to be universal at 


all latitudes. 
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Figure 2: The PD— TB, scatter plots for the squall line case shown in Fig. | from 89 GHz (a) and 166 GHz (b). The points are 
color-coded by 183.343 GHz TB. Observations over ocean are further marked by black open circles. 

The PDs observed from a conical scan are mostly positive in cloudy-sky scenes, because the cloud scattering tends to produce a 
larger TB depression due to horizontally-oriented non-spherical frozen particles. It is non-trivial to determine the magnitude of 
PD because other factors may play an important role as well. For example, as we can see from Fig. 1c, PDs are larger at 89 GHz 
over the ocean because of higher V-pol emissivity on calm water surfaces (acting like a mirror) than on windy surfaces, whereas 
land surfaces generally have little polarization due to surface roughness. In addition, liquid clouds and moisture layers may 
attenuate the PD signal from the surface (Wu and Jiang, 2004). Negative PD values can be found in the GMI data as well but 
they are mostly in the clear-sky measurements and affected by noise. The probability density function (PDF) of negative PDs at 
89 and 166 GHz behaves like a Gaussian distribution in most cases (Appendix B), from which we estimate the PD noise level to 
be ~ 0.5 K for 89 GHz and ~ 1 K for 166 GHz. It is interesting to observe that in some locations the negative PDs are non- 
Gaussian and have a significant number of measurements deviate from the Gaussian statistics. It is unclear what cause this 
statistical anomaly of the negative PDs, but it has been found in both 89 and 166 GHz data and the anomaly has an interesting 
geographical distribution. More details on the PD noise determination and the geographical distribution of the negative PDs can 
be found in Appendix B. 

The PD climatology is computed globally and sorted according to latitude, TB,, and PD value itself in 10°, 10K, and 1K bins, 
respectively. We also compiled a two-dimensional (2D) PDF with respect to PD and TB, for each latitude bin. Monthly statistics 


of the PD - TB, relationship are characterized by the peak and standard deviation of each 2D PDF. Only the cloudy-sky portion 


of the PD - TB, relationship is reported in this study using the 30 cloud detection with the 183.343 GHz radiances in each 
latitude band and computed separately for land and ocean. 

Six months of GMI L1C-R data (June, July, and August in 2014 and 2015) have been compiled to study the cloudy-sky PD 
statistics. Although we use an earlier data version (V03) in this study, we find that the results are similar between V03 and V04 
and choose to use the V03 data for the rest of paper. For the CoSSIR 640-GHz PD data, of which the campaign was based in 
Costa Rica, the statistics will be compared to the GMI data from the [0, 10°N] latitudinal bin. 

To study the scattering properties of cloud ice and precipitation-sized frozen particles at 89 and 166 GHz, we make use of 


precipitation measurements from GPM DPR in our analysis. Galligani et al. [2013] found that large 85 GHz PDs of TRMM-TMI 
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were often associated with a melting layer underneath (detected by TRMM-PR), suggesting a strong connection between cloud 
ice processes and surface precipitation. Thus, we employ the GPM Ku-band radar bright band (BB) flag to evaluate the 89 and 
166 GHz PD properties and their connection to the presence of BB, as BB is always associated with stratiform precipitation. For 
this analysis, we first collocate the measurements between the Ku-band radar and the GMI 89 and 166 GHz data, and then tag the 
BB flag to each GMI measurement. The Ka-band radar also produces a product for classifying stratiform and non-stratiform 


precipitation cases, but the Ku radar BB flag is generally considered to be more reliable (Iguchi et al., 2010). 


3. PD-TBv relationship observed by GMI and CoSSIR 


The bell-shape curve revealed in Fig. 2 is found to be universal for all latitude bins over land and ocean. Fig. 3 shows the tropical 
statistics over land and ocean for July 2014, where the PDF peaks from cloudy sky produce a similar PD bell curve for both land 
and ocean cases. The bell curve has a peak at TBy~220 K for 89 GHz and TBy~200 K for 166 GHz. Hence, despite the strong 
surface PD signal at 89 GHz the cloudy-sky PDs are quite similar between the two channels, both exhibiting a bell-curve with a 
PD peak amplitude ~ 10K. Moreover, cloud-induced brightness temperature depression exhibits a similar dynamical range 
between 100-280 K at the two frequency channels, but the same ATBy at 89 and 166 GHz features the scattering of frozen 
particles from different heights. Compared with 166 GHz, the 89 GHz TBy has a narrower dynamic range (~ 60K) before 
reaching the PD peak for both oceanic and land conditions. This result is consistent with the TMI 85 GHz PD behavior reported 
by Prigent et al. [2005], but Defer et al. [2014] reported a somewhat smaller PD peak amplitude (~ 8K) at 89 GHz using the 
MADRAS instrument on board the Indo-French Megha-Tropiques satellite mission. At 166 GHz, PD does not reach its peak 
until about 80 K colder than the clear-sky radiance. Besides, the peak amplitude at 166 GHz is slightly larger than that at 89 
GHz. This finding is consistent with Defer et al. [2014] using MADRAS’s 157 GHz polarization channel pairs. The PD bell 
curve gradually diminishes at clear-sky scenes (TBy ~ 280 K), but approaches ~ +2K at very thick convective clouds (TBy~ 100 
K). These PD-TBy features are also evident for synoptic scale clouds in the mid-latitude (see Fig. C1 or Fig. 4 for example). 

The spatial distribution of cloud occurring frequency (CF) weighted PD (i.e., CF*averaged PD in each grid box) is high over 
active convective regions (not shown), similar to the CF from tropical deep convection zones and peaking near the equator in the 
summer hemisphere. This distribution shows typical active convective centers (e.g., Central United States in July), as expected 
for thick anvil outflows from the deep convective cells. The PD-TB, bell-curve is seen at latitude bins up to 50°N/S (Fig. 4 and 
Fig. C1). 
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Figure 3: Joint PDFs of ATB - TBy statistics for the [0°, 10°N] latitude bin. The GMI 89 GHz (top) and 166 GHz (bottom) data 
in July 2014-2015 are used in the analysis, but separately for ocean (left) and land (right). Clear-sky data are excluded in 
computing the PDF peaks while the cloudy-sky PDF peaks are connected by black line with the standard deviation of PD shown 
for each bin. Cloudy-sky is determined by the “30 threshold” method using the 18343 GHz radiances. Color contour for the 
PDFs is in log scale. 

Fig. 4 summarizes the latitudinal dependence of the bell-curves at each 20° latitude bin over land and ocean respectively. Six 
months of data (June, July and August of 2014 and 2015) are averaged together to increase the statistical robustness. The surface 
PD branch has been removed using the 183 GHz radiances and the “30 threshold” separation for clear-cloudy skies. All the bell 
curves over ocean are essentially same, especially at 30°S-50°N. The wintertime high-latitude (70°S-50°S) belt is an exception. 
This is also evident at 89 and 166 GHz over land (Figs. 4c, d), which may be associated with the drier atmosphere at high- 
latitude in the winter atmosphere. Although the dynamic range of TBy shrinks with increasing latitude, the PD peaks are still 
detectable and occur at about the same TBy value (~220K) for the 89 GHz oceanic clouds credited to the large clear-sky 
variability and the resulting 30 threshold. These findings suggest that the ensemble characteristics of the ice crystal shape and 
orientation that 89 GHz is sensitive to are consistent across all latitudes, which is a piece of good news to both the satellite 
retrieval and cloud modeling group, as many cloud microphysics assumptions during the retrieval and modeling processes are 
based on mid-latitude observations. The stringency of shape and orientation parameters ensures less uncertainty in the 


retrieving/modeling of the tropical and high-latitude atmosphere. 
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Figure 4: the PDF peak curves from different latitude belts (see label) derived from June-July-August 2014 and 2015 
climatologies. Top panels for ocean and bottom panels for land conditions, respectively. Surface signal is removed using the 


183.343 GHz radiance and the “30 threshold” method. 
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The oceanic 166 GHz PDs at TBy = 200-220 K appear to increase with latitude in both hemispheres (Fig. 4b). There is also 
a hint in Fig. 4a for 89 GHz although the tendency is rather weak. This latitudinal dependence of the PD amplitude could be due 
to latitude-dependent water vapor abundance, because ample water vapor can damp the PD signal. The water vapor damping is 
stronger in the tropics and consequently flattens the PD-TB curve more. The water vapor attenuation effect is also evident for 
ocean conditions at mid-high latitude of the winter atmosphere (30°S - 70°S), which is caused by the fact that the atmosphere 
there is very dry and cannot block the surface signal anymore (e.g., Fig. Cld). However, the upper troposphere water vapor 
amount should be roughly symmetric between the two hemispheres (Solden and Lanzante, 1996), but the PD curves do not 
suggest such an inter-hemisphere symmetry. Over ocean, the PD peak amplitude is slightly larger in the summer hemisphere, 
while is contrary over land. Therefore, the role of water vapor in attenuating the PD is perhaps more complicated than expected. 


We will explore this in more detail in the next section. 
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Figure 5: Scatter plot of quality-controlled CoSSIR PD - TBy 640 GHz measurements from the TC4 campaign. Observations 
over land and ocean are not separated here. The black solid line with errorbars is the PDF peak, where bins beyond 260 K is 
excluded because of the sparseness of the data. 

Studying the CoSSIR 640 GHz measurements from the TC4 campaign, we also find a bell curve in the PD-TB relationship (Fig. 
5). Interestingly, the largest PD amplitude remains at ~ 10 K, and occurs at TBy ~ 200 K. Some high PDs are found in thin cirrus 
cloud/clear-sky scenes near TBy = 250 K and in optically thick cloud of TBy = 150K, which are also associated with large 
negative PD values. These cases are found in the July 19" and August 6" flight legs but not in July 17" flight leg. Data qualities 
are considered much noisier in the former two flights than the latter one, but we still keep to show the original data from all three 
flight legs in Fig. 5 as the peak of PD-TB relationship alters little by including the noisier data (Frank Evans, personal 
communication). 

To summarize, robust and significant PD signals from GMI’s 89 and 166 GHz channels contain ample information about ice 
microphysics of frozen particles. There exists a universal bell-curve in the PD - TB relationships at all high-frequency 


microwave channels, including the 640 GHz observations from CoSSIR. They all show a peak amplitude of ~ 10 K at TBy = 
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200-220 K. What cause(s) the bell curve is an interesting question for remote sensing microphysical properties (size and shape) 


of frozen particles. In the next section, we will discuss two plausible mechanisms that may explain the observed features. 


4. Discussions 
4.1 Radiative Transfer Models (RTMs) and the Concept of Aspect Ratio (AR) 


To evaluate impacts of frozen particle orientation and ice water content (IWC) on the TB and PD - TB, relationship, two 
numerical RTMs are employed in this study. The first model is called the Polarized Radiative Transfer Model Distribution 
(PolRadTran) developed by Evans and Stephens (1995). Two options of the scattering solver are included in PolRadTran: RT3 
and RT4. While RT3 assumes random orientation, RT4 assumes 100% horizontal alignment. The solutions are essentially 
identical for spheroid, but not for non-spherical particles. In our study, PolRadTran is configured to mimic GMI and CoSSIR’s 
viewing geometry and channel frequencies, and both RT4 and RT3 solvers are employed with assuming several non-spherical 
ice habits (snowflake dendrites, hexagonal plate, cylindrical columns, bullet-rosettes with 7 branches). Their scattering properties 
for diameter ranging from 20 to 2000 um are adapted from Yang et al. [2013], albeit the fact that Yang et al. [2013]’s calculation 
is only valid for the visible to the far-infrared spectrum strictly speaking. As described in Evans and Stephens (1995), the Fast 
Fourier Transform (FFT) approach has been used to solve the discrete dipole approximation (DDA) scattering computation 
rapidly once the ice refractive index is provided. This index has weak dependence on the ice particle temperature and we used 
the value at 240K in this study. Values at other temperatures (220K, 260K) have also been tested, and only show minor 
differences in the results. Frozen particles are assumed to obey a Gamma size distribution in general (Evans et al., 2012), which 
is characterized by mean mass-weighted equivalent sphere diameter (De). The desired gamma size distribution is achieved using 
an optimization procedure to adjust the equivalent volume diameter (D,), where Dme = (D?)/(D2). This optimization procedure 
assures that the truncated and discrete size distribution has the desired moments (i.e. Di. and dispersion width). 

The RT4 simulations we carried in this study assume a uniform ice cloud layer with horizontally oriented non-spherical ice 
particles in a background atmosphere. Lambertian and Fresnel scatterings are used for land and ocean surfaces, respectively. 
Input variables are size distribution (Dme), channel frequency, ice water path (IWP) and cloud layer bottom and top heights. 
Then, the model computes the V- and H-pol TB at the top of atmosphere (TOA) for a set of viewing angles, from which the 
GMI/CoSSIR TBs are derived by interpolating the model TB outputs to their viewing angles. 

The RT3 option of PolRadTran model will also be employed to reveal the radiation behavior of a medium of randomly 
oriented particles. Different from RT4, which only considers the thermal radiation, RT3 includes the solar radiation as well. But 
the for our purpose on the microwave spectrum, solar radiation is negligible. 

The RT4 model and the configuration of our simulations have several limitations. Firstly, PolRadTran is not a fully-polarized 
RTM, strictly speaking. Plane-parallel, horizontally homogeneous and azimuthally symmetric orientation are assumed for 
PolRadTran and RT4, such that we only need to calculate two out of the four Stokes parameters while the other two remain zero. 
The modeled ice crystal shapes (except dendrites) represent primarily the pristine mid-latitude cirrus from synoptic outflows, 
which is not ideal for larger ice particles in the stratiform precipitating regions nor for tropical deep convection cases. The ice 
particle scattering properties were adapted from Yang et al. [2013], where the calculations stop at far-infrared. In the future, more 
suitable database, e.g., Liu (2008), might be employed to replace the current input database. While the DDA calculation is very 
efficient for computing cirrus particle scattering, it is still computationally expensive for large particles, of which the valid range 
of IWP is limited to <1 kg/m? and De to <500 um. Besides, the RT4 surface model and homogeneous ice cloud layers are often 


over-simplified compared to the reality. Because of these limitations, RT4 is perhaps mostly suitable for simulating the 640 GHz 
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PDs where the surface contributes the least, compared to the lower frequency channels, and this channel is insensitive to 
precipitation-sized ice particles. Assuming randomly particle orientation, Evans et al. [2012] applied the RT4 model with a 
dynamical mixture of different particle shapes to cloud ice retrievals on the CoSSIR data. 

Another RTM for simulating frozen particle scattering, called Cloud Radiance Model (CRM), is based on the cloudy-sky 
forward model in Aura Microwave Limb Sounder (MLS) (Wu and Jiang, 2004; Gong and Wu, 2014). CRM is a non-polarized 
model, and it assumes only spherical ice particles in the scattering calculation. It has been applied to simulate cloud-induced 
radiances at frequencies < 3 THz for a complicated viewing geometry (e.g., satellite nadir/limb sounding, in-cloud aircraft flight, 
and ground stations). Like RT4, CRM assumes the Gamma size distribution in the current study. Different from RT4, CRM 
allows the user to specify the vertical distribution of IWC or liquid water content (LWC) in the mixed-phase case, to study 
sensitivity of the cloud-induced radiance to cloud layer height. Cloud height is another free parameter in addition to particle size 
and IWP, which can alter the cloud radiance amplitude significantly in similar frequency range (Gong and Wu, 2014). 

To study the effect of frozen particle orientation with CRM, we assume slightly different profiles of scattering volume 


coefficients at V- and H-pol, ty(z) and t,(z), respectively). The difference is characterized by the aspect ratio (AR) factor, 


defined as AR = “4 (2) )> which is assumed to be independent of height. As one can see from the definition, AR is a 


Ty(z 
function of all microphysical property that plays a role in determining the optical depth along the line-of-sight (LOS), including 
particle size, orientation, habit, density, etc. Only by assuming all other ice microphysical parameters are homogeneous along the 
LOS can we claim that AR directly corresponds to the axial ratio of the major and minor axes of the ice particle projected to the 
LOS. The definition of AR is analogous to the “dichroism factor”, but geo-magnetic field is believed to play the deterministic 
role for the latter (Mishchenko, 1991; Mishchenko et al., 2002; Yang et al., 2011). In simulating the PD-TB relationship, we vary 
the AR value but keep the rest model input parameters (e.g., Dime, IWC profile, etc.) unchanged. This is equivalent to the particle 
AR effect in which horizontally-oriented particles tend to create a stronger scattering for the H-pol radiation than for the V-pol. 
Thus, we carry out two parallel simulations with different cloud scattering coefficients but identical settings for other parameters. 
We then associate the two simulations with the V- and H-pol cases to evaluate their differences or PD as would be observed by 
satellite/aircraft at a slant-viewing angle. Using this CRM approach, we will also be able to study the environmental effect on the 
expected PD as a function of water vapor attenuation in the atmosphere. 

In the next section, we will proceed the model explanation from an extremely simplified two-layer model. By computing the 
layer-by-layer radiative transfer with including the AR concept, we can reproduce the bell-curve with reasonable range of PD 
values. Then, the more sophisticated RTMs described above will be employed for further simulating and understanding the 


observed PD - TB characteristics. 


4.2 A Simple Model for the PD-TB relationship 

The leading hypothesis of the observed positive PD in cloudy sky is that the bulk ice properties produce different volume 
scattering coefficients in the V- and H-polarization. The bulk volume scattering coefficients can differ between the V- and H- 
polarization, depending on ice crystal shapes and their orientation. While the bulk scattering effects of cloud particles are 
detectable, details about how particles of various shapes organize themselves inside a cloud are beyond the capability of current 
remote sensing techniques. Clouds may produce little or no PD if the shaped particles are randomly oriented (e.g., tumbling in a 
turbulent environment). On the other hand, as discussed below, the PD tends to decrease with cloud optical thickness as a result 
of cloud self-extinction in the radiative transfer. 

Let us assume that frozen particles are horizontally aligned, and the V- and H-pol extinction coefficients are dominated by 


scattering. Naturally V- and H-pol radiances pass through the same cloud with slightly different optical depths (ct, and ty). We 
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define the aspect ratio (AR) factor as AR = t,/Ty to characterize this optical depth difference. Fig. 6 illustrates a simple two- 
layer model, where cloud optical depth is 77. Neglecting the water vapor attenuation, radiance at the top-of-atmosphere (TOA) 
can be expressed analytically as, 

TBy =Tyye °Y +T,(1—-e-?Y) (1) 

TBy = Tyye 4 + T,(1—e 2) (2) 
where T, is the cloud source function, T,y and T,,, are the polarized radiances from the background below the cloud resulted 
from the combination effect of surface, rain, or liquid cloud emissions. The cloud source function T;= @9 *Tsea: + (1-@0)'T> 
contains cloud emission (/-w,)-T, and scattering @°Tsca; contributions, where wa is single scattering albedo, and T> is the cloud 
emission at local air temperature. T,.q, is the radiance scattered by clouds into the field of view, a result of the scattering phase 
function convolved with surrounding radiation, which is approximately half of the upwelling radiation if the downwelling 
radiation is very small (Wu and Jiang, 2004). The factor (1 — e~*2) takes into account the self-extinction process that is intrinsic 
to any source function in the radiative transfer. In the case that the background emission is unpolarized, i.e., Tyy = Tyy = T1, 


then from Eqs. 1-2 we have 


PD= TBy _ TBy = (T; = T,;)(e7 "2 = e *2H) (3) 
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Figure 6: Schematic two-layer model with a liquid/rain/surface layer at T,; and ice cloud layer at T2 to explain the observed bell 


curve in the PD-TB relationship. See the text for more details and definition of the variables. 


This simple 2-layer model (Eqn. 3) is able to reproduce several fundamental features of the observed PD-TB relationship. 

As shown in Fig.7a, it predicts the bell curve as seen in the GMI and CoSSIR polarimetric measurements, showing small PDs at 

both small and large cloud optical depths or respectively warm and cold TBs. At large optical depths, as implied by Eqn. 3, the 
PD diminishes. At small optical depths, Eqn. 3 can be approximately written as 

PD = (T, — T) (2H — Tay) (4) 

Because T; — T; is roughly constant according to the observation shown in Fig. 3 in the tropics, this linear dependence on T27, — 


Toy helps to spread the slopes for different ARs because T24 — Tzy = (AR — 1)t2y. The model can also explain the complex PD- 
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TB relationships seen in the GMI oceanic data where the background surface emission is polarized. Because the 89-GHz channel 
is more sensitive to the polarized surface emission than the 166-GHz, the observed cloud PD often contains both surface and 
frozen particle contributions. As shown in Fig. 7b, the combined PD from the surface and cloud scattering can result in a tilted 
bell curve. The observed PD statistics in Fig. 3 are likely a mixture of the two cases shown in Fig. 7. 

Finally, the conceptual model predicts weak dependence on channel frequency for the PD-TB relationship. The reason is as 
follows. If cloud extinction is dominated by scattering, ct is proportional to the 4" order of the frequency, ie., T= af*. 
Therefore, AR is independent of frequency by definition. If the background emission is unpolarized (i.e., Tjy = T14 = T,), the 
PD from Eqns. 1-2 has a maximum of (T, — T2)(e~*2” — e~4®*"2V) at t2y = In (AR)/(AR — 1), which is only weakly 
dependent on frequency. This weak-dependence of PD on frequency seems to be consistent with the observations that PD 


amplitudes barely increase from 89, 166 to 640 GHz. 
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Fig.7. PD -T, relationships from Eqns. 1-2 for different aspect ratios (ARs, value marked above each curve), assuming that the 
background emission is unpolarized (left, Tjy - Tjy=0) and polarized (right, Tjy - Tjy=10K). The dashed line indicates a 
typically polarized surface emission from wind-induced surface roughness. 

Radiative transfer in this conceptual model is highly simplified. Although it greatly facilitates our understanding of the 
bell-curve feature of the PD - TB relationship, it does not necessarily yield any information about ice microphysical parameters, 
and hence is not very useful for the development of retrieval algorithms. In order to test if this theory works in more realistic 
situations, and more importantly, in order to identify parameters that contribute and are sensitive to the PD — TB relationship to 
shed light on future retrieval algorithm development, we further test the AR values of 1.2, 1.3 and 1.4 using the non-polarized 
CRM. This model is detailed in Wu and Jiang (2004). To simulate polarized radiances, the CRM is run twice: one with the 
regular cloud extinction profile and one with the regular profile multiplied by AR. The two simulated radiances are differenced 
to yield the PD as a function of AR. The surface temperature is set to 300 K in all simulations. The vertical distribution of IWC 
is assumed as a typical deep convective cloud profile with an anvil top (see Wu and Jiang, 2004 for a description), where the 
cloud bottom is 5 km (right above the freezing level in the tropics) and tops at 13 km. Surface polarization is neglected in the 
simulations. A standard tropical atmosphere is used as the clear-sky background, and hence effects of the water vapor attenuation 
are included. Gamma distribution is used for ice particles with De = 200 pm. This large value takes into account a mixture of 
cloud and precipitation particles such that large TB depression can be produced. Other Di. values have also been tested. While 
the dynamic range of TB is sensitive to Die, the peak PD amplitude is not. In the simulations only IWP is varied to create the PD- 


TB curves that best mimic the observations. 
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Figure 8: CRM simulated PD-TB curves in black for AR values of 1.2, 1.4 (thin) and 1.3 (bold), overlaid on top of the GMI and 


CoSSIR observations in red for (a) 89, (b) 166 and (c) 640 GHz with the standard deviation shown in pink shades. Surface 


contaminations are removed using the cloud flag derived from the 183.343 GHz channel. 

Although AR in the range of [1.2, 1.4] produce a good envelope for the observed bell curves, AR = 1.3 seems to fit best overall 
to the observations. This AR value was also found in Davis et al. (2005) in studying the 122 GHz data, but AR defined in that 
paper has a direct microphysical meaning of particle axial ratio. Having said that, there is an indication that the best-fit of AR 
might increase with frequency, suggesting an increasing contribution of horizontally-oriented particles at higher frequencies. 
Since the higher frequencies are more sensitive to small frozen particles, it implies that small frozen particles have an increased 
percentage of the horizontal orientation. While it is difficult to infer the detailed composition of particle habits in these clouds, 
ice crystals in cirrus clouds (more sensitive by 640 GHz) are likely to have a predominant type (e.g., plate, column, etc., 
Libbrecht, 2005) and snow aggregates (more sensitive by 89 GHz) tend to have diverse irregular shapes. In addition, the ambient 
air in the upper troposphere anvils and cirrus (where 640 GHz weighting function peaks) is likely more stable than the tumbling 
environment inside deep convective cores (where 89 GHz weighting function peaks). Therefore, higher degree of alignment of 
non-spherical particles in the anvils/cirrus could also contribute to the larger AR value at 640 GHz. This possibility will be 
elaborated further in the end of Section 4.3. The third possible explanation involves with the water vapor attenuation effect. 
Although all three frequencies are window channels, they are still affected by water vapor abundance. But as found previously in 
Gong and Wu (2014), the water vapor attenuation effect is minor when cloud is present. The last plausible explanation of why 
the simulated 640 GHz peak amplitude is weaker than that of the other two frequencies with the same AR value originates from 
the background atmosphere lapse rate change. As the atmosphere becomes more isothermal approaching to the tropopause, the 
differentiation between the slightly higher weighting function (WF) of the V-pol compared to that of the H-pol starts to diminish 
accordingly, resulting smaller PD at 640 GHz as its WF peak centers at the upper-troposphere. Although the simulated behavior 
of the peak amplitude change with frequency is likely due to the above reason, the observed peak amplitude at 640 GHz is not 
discernably weaker than the other two frequencies. As one can always argue that this might be due to limited sample size for the 


640 GHz measurements, including more 640 GHz observations in the analysis indeed is the only way to solve the puzzle. 


4.3 Horizontally Polarized Model Simulation 


AR is a parameter that characterizes the collective scattering property of frozen particle habits. To study the relative 
importance of individual particle habit, we use a horizontally polarized RT4, a variant of the PolTranRad model [Evans and 


Stephens, 1995], to simulate the PD-TB relationship. Four types of ice habits and different D,,. are used in the simulations, and 
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the results are shown in Fig. 9 for 640 GHz assuming all frozen particles are horizontally aligned. As seen in Fig. 9, snowflake 
dendrites and cylindrical columns are most effective shapes to produce large PD values at 640 GHz, whereas cylindrical columns 
and bullet-rosettes can also induce a similar PD amplitude to what is measured by CoSSIR. The first two ice crystal types result 
in excessively large PDs in the PD-TB relationship. For bullet-rosettes, on the other hand, the PD does not produce the bell curve 
when the cloud becomes optically thick, whereas the cylindrical column shaped crystal with Dre between 200 and 300 um 
produces the best agreement with the observation. According to Yang et al. [2013], this type of ice crystal with Dine = 250 um, 
when projected to GMI’s viewing direction, has an AR of 1.46, a value in a surprising good agreement with that from the CRM 
simulations shown in Fig. 8c. This consistency also supports the AR values used in the simple conceptual model and its radiative 
transfer argument for the bell-shape curve. We hope with more accurate multi-frequency polarimetric measurements that ice 
habits may be retrievable from the PD - TB relationship eventually. More questions also rise up. For example, given that ice 
crystals have a vast variety of shapes in reality, why are the PD - TB relationship at 640 GHz similar to those at lower-frequency 
channels? Future field campaigns with polarimetric measurements from higher frequency channels should provide more 


insightful information to address similar questions. 
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Figure 9: PolTranRad simulations of 640 GHz PD - TB relationship for four different ice habits: (a) snowflake dendrite, (b) 
hexagonal plate, (c) cylindrical column, and (d) bullet-rosette with 7 branches. All frozen particles are assumed to be 
horizontally aligned. The red and blue in the snowflake dendrite simulation correspond to Dye = 200 and 300 um, respectively. 
In the rest panels, red, yellow, cyan, and blue denote Dye = 200, 300, 400 and 500 tum, respectively. 

One of the top questions investigated in this study is: what process causes the low PDs of the bell curve at cold TBs? One 
would argue that at cold TBs (i.e., strongly depressed by the scattering of deep convective clouds) the orientation of frozen 
particles were mostly random instead of horizontally aligned. As a result, those deep convective cold TB cases would yield a 
weaker PD. PolTranRad has another version called RT3, which allows to simulate effects from randomly orientated ice crystals, 
and the bell curve is found to be still present, despite the amplitude is rather small (< 3K) regardless of type of ice habit used. 
Thus, the random orientation of frozen particles does not completely cancel out the particle shape effects, as revealed by the RT3 
simulations. 

A varying mixture of randomly and horizontally-oriented particles may also provide a plausible explanation of the observed 
bell curve. If the larger mixing factor “a” corresponds to a more mixed environment inside clouds, o is likely to increase with 


deep convection strength, or a colder TB. On the other hand, horizontally aligned non-spherical particles are likely more 
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effective to produce PDs in the anvil outflows than convective cores, because random orientation within the turbulent convective 
core tends to cancel out the PD. This mechanism seems to have a support evidence in the ground-based Doppler radar 
observations (Homeyer and Kumjian, 2015), where they often found lower differential reflectivity (ZDR) near the deep 
convective core but larger values in the anvil outflows. In addition, the irregular shape of graupel found near convective cores 


likely contributes to a small PD as well as suggested in their paper. 


4.4 PDs from melting layers and ice clouds 


Since 89 and 166 GHz radiances are sensitive to scattering and emission of the precipitation-sized frozen particles, we need 
also to take into account the effects of the floating snow layer (Xie et al., 2012) or the melting layer close to the surface 
(Galligani et al., 2013) underneath the cloud layer. Although it is difficult to quantify the PD sensitivity separately for different 
layers of frozen particles, it is possible to evaluate the response of PD to the melting layer with the combined DPR-GMI 
observations. Melting layers manifest themselves as a unique bright band (BB) on the Ku-band radar reflectivity profile, as a 
result of stratiform rainfall. We therefore use the BB flag and precipitation flag in DPR’s Ku-band radar Level-2 products to 
discriminate the stratiform and non-stratiform precipitation scenes. Non-stratiform precipitation types include convective rainfall 
and snowfall. 


(a) Ocean, PD [kK] (b) Land, PD [K] 
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Figure 10: Joint PDF of 89 and 166 GHz PDs during January 2015 over (a) ocean and (b) land for precipitation with the bright 
band (color shaded) and without the bright band (black contours). Color and contours are normalized such that the total areas 
equal to unity, and are plotted in log scale. 

Fig. 10a shows a joint 89 and 166 GHz PD PDF with rather diverse relationships for oceanic measurements with and 
without BB. There are two branches in the 89 and 166 GHz PD correlation characterized by different slopes: 1:2 and 4:1. The 
measurements in the 4:1 slope are dominated by the polarized ocean surface emission, whereas the 1:2 group has little 
contributions from the surface. Because the 1:2 branch over ocean is also associated with the strongly polarized ocean surface, it 
should be related to lightly precipitating scenes that otherwise would show much stronger PD signals. In the case with BB, most 
of the polarized surface emission is blocked by BB. For non-stratiform precipitation (without BB), most likely there will be no 


PD signal from both channels (i.e., the center of the black contours is near the origin point in Fig. 10a). This indicates that the 
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dominant PD signals of 89 and 166 GHz indeed come from ice cloud for non-stratiform precipitation, even it is snowfall when 
the snow scattering is relatively weak at 89 GHz, which makes it capable of “seeing” the surface in this case. When the melting 
layer is present, both 89 and 166 GHz PDs incline to be positive. As a result, when the melting layer is present, the peak of the 
joint PDF (dark red) moves to ~ 5K and 8K for 89 and 166 GHz, respectively, as shown by the color shades in Fig. 10a, implying 
that the melting layer in fact contributes partially to the PD signal. This finding is consistent with Galligani et al. (2013) where 
they found that 85 GHz PD of TRMM’s TMI is likely positive when the melting layer is present. It is also consistent with the 
higher BB ZDR as observed by ground-based radars (e.g., Oue et al., 2014). While 166 GHz is more sensitive to smaller ice 
crystals than 89 GHz, the more positive PD response in the 166 GHz indicates that smaller ice crystals should be more obsolete 
and/or more horizontally oriented. Larger particles like snow aggregates likely exhibit irregular shapes and the scattering signal 
in 89 GHz would hence be less directional, resulting a smaller PD response in this channel. This finding is highly consistent with 
the weak increase of AR value with the increase of channel frequency in CRM simulations shown in Section 4.2. Moreover, 
frozen particles of all sizes tend to be more horizontally oriented for stratiform rainfall compared with non-stratiform 
precipitation (Houze, 2014, Section 6.2). 

Over land, the PD characteristics vary roughly the same way with that over ocean, but the 4:1 branch exhibited in the ocean 
case is not present in the land non-stratiform precipitation cases. The PDF peaks along the 1:2 line are different in the cases with 
and without BB. Other factors, e.g., liquid emission, may also play a non-trivial role here, but they are hard to be justified at this 
moment. Due to the complexity of melting layer, this paper does not intend to further quantify its role on the PD’s 


characteristics. 


5. Conclusion and Future Work 


High-frequency microwave polarimetric measurements contain ample information of frozen particle habit and orientation. In 
this paper, the GMI 89 and 166 GHz PD characteristics are analyzed, which reveal a universal bell-curve in the PD - TB 
relationship. The bell curve also exists in the CoSSIR 640 GHz observations. All three frequency channels show a maximum PD 
amplitude of ~ 10 K, occurring around TB = 200 — 220 K. Such a PD-TB relationship is robust for almost all frequencies and 
latitudes. Moreover, both 89 and 166 GHz channels show positive PD responses when a melting layer is present. This is a new 
feature not reported before, suggesting that at least part of the PD signals comes from the melting layer. The melting layer 
contribution is confirmed by comparing the PD statistics with non-stratiform precipitation situations. 

Using a simple analytical radiative transfer model, we can explain the bell curve of the PD - TB relationship with the term of 
aspect ratio (AR). In this simple model we assume that the V- and H-pol extinction coefficient profiles are similar and scaled by 
AR. According to this conceptual model, the PD is resulted from the stronger extinction at the H-pol but diminishes as the cloud 
scattering is too strong to become saturated to a constant TB in both polarizations. In addition, a fully polarized RTM is 
employed to study the effects of the frozen particle habit and orientation on the PD. The mixing of different particle habits and 
orientations inside deep convections break down the directional orientation, which cancels out the non-spherical shape effect, 
which can result in a bell curve in the PD-TB relationship. 

Our fully polarized RTM simulations suggest a strong sensitivity of the PD’s magnitude to ice habit at 640 GHz, i.e., the 
shape factor. Previous research (Defer et al., 2014) using another fully polarized RTM also recognized that other ice 
microphysical properties, including density, size, etc., would impact the simulation results seriously. Therefore, more 


observations from higher-frequency microwave or even IR channels are very much needed, and would greatly facilitate the 
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decode of the microphysical information of frozen particle shape and orientation at different portions of the crystal size 
spectrum. 

Negative PD values are also found in these polarimetric measurements. The spatial distribution of the negative PDs does not 
coincide with the most active convection regions. Although they have been hypothesized as a result of vertically-oriented frozen 
particles, careful consideration of measurement noise and surface properties is required before interpreted as a result of particle 
properties. The ample information embedded in PD sheds light on the future path of retrieval method development of more 
frozen particle microphysical properties. 

Last but not the least, the observed PD-TB relationship has an important implication for cloud ice retrieval. Gong and Wu (2014) 
used an empirical IWP - TB relationship derived from CloudSat-MHS (Microwave Humidity Sounder) measurements for the 
IWP retrieval, where they found this relationship is nearly linear for medium thick ice cloud (i.e., anvils). The observed PD value 
range in this study therefore can be translated into a 30% IWP retrieval error if polarization is neglected. This is a very rough 


estimate that warrants a thorough evaluation in the future. 


Acknowledgments 


We thank Frank Evans for providing the RT4 code with GMI configuration, and help on interpreting the CoSSIR data. We are 
also grateful to Ben Johnson, Stephen Munchak and Paul Racette for helpful discussions. Jeffrey Piepmeier and Rachael 
Kroodsma are acknowledged for assisting the interpretation of GMI’s channel noises. Two anonymous reviewers’ and Yana 
Mentrok’s insightful comments and suggestions are highly appreciated. This work is supported by NASA NNH12ZDA001N- 
INVEST fund. 


19 


10 


13 


20 


Appendices 


Appendix A: Mismatch and re-collocation of GMI low- and high-frequency channel measurements 
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Figure Al: Scatter plots of 1-day tropical oceanic radiance measurements from L1B (a and c) and LIC-R (b and d) datasets, 


showing poorly (L1B) and correctly (L1C-R) georegistered radiance data. Points are color-coded by the 183.343 GHz radiance 
values with red being warm and blue being cold. 


Fig. Al illustrates the necessity of using L1C-R data to screen out the 89 and 166 GHz clear-sky scenes. Compared to LIB 


dataset that does not collocate the LF and HF channels, 183.343 GHz TB now corresponds much better with 89 GHz TB. 
Moreover, 166 and 89 GHz TB shows a much better coherence in the L1C-R data, showing two separate branches more clearly 


in Fig. Al (d). The different slopes of these two branches correspond to the clear-sky and cloudy-sky characteristics, 


respectively, and the clear-sky branch also is the warmest in 183.343 GHz TB. Scatter behavior of Fig. Al(d) is very similar to 


comparable instrument, e.g., Special Sensor Microwave Imager and Sounder’s 91 and 150 GHz channels (not shown). 


Appendix B: Estimation of Instrument Noise and Distribution of Negative PD 


GMI’s 89 and 166 GHz PDs are not always positive. As a matter of fact, negative values occur frequently, and their 
occurrence is independent with TB. However, large negative PD only happens occasionally. The histogram of negative PD 
closely follows the Gaussian distribution, as shown in Fig. B1 by the black solid fitted lines. The “Gaussian_fit” function from 
IDL software package is used for this purpose, where the mean (1) and standard deviation (o) parameters are given as output. 
The standard deviation of the fitted Gaussian curve is a good estimate of the instrument noise (e.g., Wu et al. [2009]). Following 
the same analogy, we define the standard deviation parameter from the fitted Gaussian distribution as the GMI’s instrument 


noise. G is computed consequently for two months of 2015, January and July, for each 10° latitude belt, which varies limitedly 
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among different latitudes, surface condition and season. The averaged o is about 0.5 K for 89 GHz and | K for 166 GHz. Note 


that the radiance noise level should be derived by dividing o by V2 . 

If the negative PDs are purely originated from instrument noise, they should be randomly distributed geographically. 
However, that is not the case for 166 GHz. If we plot out the occurrence frequency (OF) map of negative PDs that are smaller 
than the negative instrument noise estimated using the aforementioned method, the OF is trivial at 89 GHz, indicating that 
negative PD is likely instrument noise. But for 166 GHz, large negative PDs occur quite frequently (up to 18% of the time), and 
they seem to favor the winter-hemisphere side of the tropical deep convections, as shown by Fig. B2. This is the first time such a 
geographical preference has ever been reported. 

One theory explains the negative PD as the scattering effect from predominantly vertically oriented ice crystals that may be 
oriented by the electromagnetic field created by lightning near the cloud top (Prigent et al., 2005; Homeyer and Kumjian, 2015). 
We can also consider the vertical orientation a counterpart solution of the horizontal orientation from the RTMs as the latter can 
only produce positive PD signals. However, global lightning distribution, which essentially connects to the tropical deep 
convective activities, correlates poorly with the geographic distribution of negative PDs (e.g., Cecil et al., 2014). Instrument 
noise is not likely the underlying cause for 166 GHz either since it would be randomly distributed otherwise. Having said that, 
the secondary branch in the negative PD branch in Fig. 4 surely raised up more interests on studying the causes of negative PD. 
Although the actual cause remains unclear at this moment, negative PD, especially the large magnitude one, assures a rather 


attractive topic for future studies. 
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Figure B1: Histogram (cross) of negative PDs based on all GMI observations during January 2015 within the latitude belt of [0°, 
10°]N for 89 (top) and 166 GHz (bottom) for ocean (left) and land (right) conditions. The Gaussian distribution fitted curve is 


plotted in black solid line, and the standard deviation from the fitted curve is marked on each plot. 
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Figure B2: Geographical distribution (colored contours) of negative PD occurrence frequency during January (a) and July (b) 
2015 for 166 GHz. Gray shades marked the occurrence frequency of TB, < 160 K scenes, which is essentially equivalent to the 
occurrence frequency of deep convections, except that the signal over Tibet Plateau in Fig. B2(a) is probably induced by surface 


contamination. 


Appendix C: PD — TB relationship at high latitudes 


Fig. Cl shows the 2D PDF distribution from the same month climatology with Fig. 3, except that it is from another mid- 
latitude bin. Compared with Fig. 3, one can see very similar behavior of the PDF peak-connected lines. But the surface signal 
starts to emerge at 166 GHz at mid-latitude, as shown by the warm branch in the bottom panels of Fig. Cl. “30 threshold” 
method is proven here to be able to effectively exclude the surface signal in the PD — TB relationship, although it might be too 


stringent, especially for 166 GHz, as the clear-sky PD is only weakly dependent on the ocean surface roughness for this channel. 
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Figure C1: Same with Fig.3, except for [45°, 55°N] latitudinal belt. 


Appendix D: Cloud Occurring Frequency (CF) derived from the “30 method” 


To demonstrate the credibility of the “30 method” on screening out the clear-sky scenes, we show the CF of July 2015 from the 
89 and 166 GHz channels in Fig. D1. The geographic distribution of CF of anvil and deep convective clouds agrees well with 
that observed by the CloudSat radar (Sassen et al., 2009). Due to the stringent criteria set by the “30 method”, many cirrus scenes 


are classified as clear-sky scenes, and therefore cirrus cloud CF is obviously biased low using this method. 
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Figure D1: Cloudy-sky cloud occurring frequency (CF) by applying the “30 method” from GMI’s 89 and 166 GHz observations 


during July 2015. 
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